Multiple - vr transitions in SIFS Josephson tunnel junctions 
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We report on experimental studies of superconducting coupling through a thin Ni76Al24 film. A 
new patterning process has been developed, which allows in combination with the wedge shaped 
deposition technique the in situ deposition of 20 single Nb/Al/Al203/Ni3Al/Nb multilayers, each 
with its own well defined NiaAl thickness. Every single multilayer consists of 10 different sized 
Josephson junctions, showing a high reproducibility and scaling with its junction area. Up to six 
damped oscillations of the critical current density against F-layer thickness were observed, revealing 
three single 0-7r-transitions in the ground state of Josephson junctions. Contrary to former exper- 
imental studies, the exponential decay length is one magnitude larger than the oscillation period 
defining decay length. The theoretical predictions based on linearized Eilenberger equations result 
in excellent agreement of theory and experimental results. 

PACS numbers: 74.25.Ha, 74.50.+r, 01.30.Rr 



In the past few years there was a noticeable interest 
to the unconventional Josephson junctions 0, Q, 01; in 
particular, to the so-called tt - junctions having negative 
critical current. These junctions provide a tt - shift in 
the ground state and were realised experimentally in SFS 
(supcrconductor-ferroniagnet-superconductor) and some 
HTS structures. 

The intensive experimental study of "0" -"tt" transi- 
tion in SFS Josephson junctions !1 i, 1 S S i, ES [H 
Ha . Il3j | confirms the existence of critical current oscilla- 
tions upon the thickness of ferromagnetic interlayer df. 
Different structure of SFS sandwiches and SIFS tunnel 
junctions having been fabricated up to now. They con- 
tain regions which are controlling the critical current and 
difficult to control in experiment and describe in theory. 
They are SF interfaces, dead layers, and the region in 
S banks with suppressed superconductivity. Contrary to 
that the bulk properties of F material can be better con- 
trolled and well described by theoretical models based 
on quasiclassical theory of superconductivity. These the- 
ories predict that for large thickness of ferromagnet the 
critical current of SFS junctions have to exhibit a damped 
oscillations as a function of df 



effective lengths. They are the decay length ^i, and the 
length ^2, which determines the period of oscillations. In 
dirty limit the expressions for ^12 follow from the Usadel 
equations |l4| and have the form 
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where Df and Eex are the diffusive coefficient and ex- 
change field of ferromagnetic material, respectively. In 
the clean limit one can easily get from Eilenberger equa- 
tions that 15] 
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Here di is the position of the first minima, Iddo) is the 
first experimental value of Ic{df). These two values take 
into account the resultant action of SF interfaces and 
their vicinities. The oscillations are characterized by two 



where I ^ ^0 is the electron mean free path and vp is 
the Fermi velocity in a ferromagnet. 

It is clearly seen from Eqs. and Q that for dirty 
materials ^2 > S.i, and in the limit of large exchange 
energy, Eex >> ttT, the characteristic lengths are nearly 
equal ^1 — ^2- In the clean limit these ^1 and ^2 are 
completely independent. 

Our analysis of both the bulk properties of ferromag- 
net material s fla a nd the experimental data [1, a S S 
II H in in ElEiLlias shown (see Table Hj) that in di- 
lute ferromagnets 0, H, K 113 the electron mean free 
path is very small providing the fulfilment of the dirty 
limit conditions for the F interlayer. In these experi- 
ments ^1 ~ ^2 as it follows from Eq. ((SJ. Contrary to 
that, in the structures with Ni interlayer the relation be- 
tween ^2 and ^1 is just the opposite so that more complex 



TABLE I: Characteristic lengths in ferromagnetic materials 
for SFS Josephson junctions. 



Ref. 


5i(nm) 


^2(nm) 


F-material 


VF (m/s) 


Ee^K) 


[8] 
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1.6 


Fe2oNi8o 


2.2 • 10^ 


1100 


[11] 
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2- 10^ 


400 


[10] 
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3.5 


Cuo.53Nio.47 




850 


[5] 






Cuo.52Nio.48 






[9] 


1.7 


1 


Ni 


2.8 ■ 10^ 


2300 


this work 


4.6 


0.45 


NisAl 


1.5 ■ 10^ 


1000 



model [131 should be used for the data interpretation. It 
is necessary to point out that in all previous experiments 
(except [13) the structures were not fabricated in one 
run, so that the certain degree of non-reproducibility of 
magnetic constants of F materials occurred for junctions 
with different df. This results in increase of spread of 
data with increasing df and did not permit to observe 
the large amount of oscillations. 

In this work, we improved the reproducibility of the 
junction parameters by preparing the structures in one 
run with F - layer thickness between 10 nm and 20 nm. 
For the first time, we succeeded in observation of up to 
six damped oscillations of critical current with F - layer 
thickness. To do this we have used the "wedge" shaped 
F - layer technique and a new ferromagnetic material 
(NiaAl). We experimentally obtained a one order of mag- 
nitude difference between ^1 and ^2 which is consistent 
with theory based on the Eilenberger equations [13 • 

The bottom electrode of SIFS samples consists of 
Nb/Al/Nb/Al and was deposited on oxidized 2 inch Si 
wafers with argon magnetron sputtering. The top 10 nm 
thick Al layer has been oxidized for 2 minutes in a 
0,1 mbar pure oxygen. The following NiaAl interlayer was 
sputtered with neon gas from one single target. The tar- 
get composition of the alloy was determined by Ruther- 



ford Backscattering (RBS) and is Ni74Al26. A 30 nm 
Nb top layer was deposited in situ to prevent the inter- 
layer from oxidation. A schematic cross section of the 
deposited multilayer is shown in insert of Fig. 1. To 
achieve a thickness - gradient of the NisAl layer the po- 
sition of the target above the substrate has been shifted 
during deposition for several centimeter. This permits 
us to produce a rather linear thickness gradient over the 
whole 2 inch substrate (see Fig.QJ. 

The transport measurements were performed at 4,2 K. 
Figure [3 shows the dependence of the critical current, 
Ic, upon external magnetic field H. The rather optimal 
agreement between experimental results and the Fraun- 
hofer function fit indicates a uniform and homogenous 
current distribution in the junction. The current- voltage 
characteristics (CVC) of several Josephson junctions can 
be seen in Fig.|21 They show a clear superconducting tun- 
nelling behaviour with a hysteresis in its curves. These 
four curves belong to one defined NisAl layer thickness of 
12.5 nm. The current is normalized to the current den- 
sity, since the four junctions differ in their sizes. It can 
be seen in Fig. |3| that the current density differs only 
several percent, thus showing a good reproducibility of 
our junctions. This result also indicates that within one 
patterned line the F - layer thickness variation could be 
neglected. The thickness gradient has been characterized 
and proofed with RBS measurements. The variations are 
several Angstroms. We developed a patterning process 
that allows the creation of 20 different separated 500 /im 
wide lines, distributed homogeneous along the F - layer 
thickness gradient over the 2 inch wafer. Each lines con- 
sists of 10 different sized circular Josephson junctions. 
The junction area differs from 5 to 1000 /.tm^. The varia- 
tion of F - layer thickness inside an individual junction is 
negligible small. The ferromagnetic properties of NiaAl 
films depend on neon pressure. Details will be published 
elsewhere 18]. The magnetic properties of the NiaAl lay- 



NijAI samples 
linear fit 




10 15 20 25 30 35 40 

relative position on wafer, mm 

FIG. 1: Ni76Al24 film thickness measured at different posi- 
tions on a 2 inch Si wafer with the RBS method. Inset: layer 
sequence in a cross section. 
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FIG. 2: Fraunhofer - pattern of a typical circular 
Nb/Al/Al203/Ni3Al/Nb Josephson junction, measured at 
4.2K (rectangles). The junction size is about lOOO/im^. The 
solid line results from a fit to the Fraunhofer function. 



600 
400 
200 

^ 

o 

-200 
^00 
-600 



■-1000|jm' 

■ - 750 \ini 

■ - 500 [jm" 
■-lOOMm' 




-0,5 -0,4 -0,3 -0,2 



-0,1 0,0 0,1 

U, mV 



0,2 0,3 0,4 0,5 



FIG. 3; Current density vs. voltage characteristics of four dif- 
ferent Josephson junctions from one patterned line with dif- 
ferent junction sizes, measured at 4.2 K. The size of junctions 
changes from 100 ^,m? to 1000 /im^. The F-layer thickness is 
12.5 nm. 
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FIG. 4: Critical current density of up to 60 single SIFS 
Josephson junctions against the F - layer thickness. The solid 
line indicates a theoretical fit after Eq. Q with the two fitting 
parameters ^i = 4.6 nm and ^2 = 0.45 nm. 



ers were measured with a SQUID - magnetometer. The 
Curie temperature for a 250 nm thick NisAl layer is 74 K. 
The Nia Al thickness dependence of the critical current 
density can be seen in Fig. 0] Several Josephson junc- 
tions of each line were measured and plotted versus their 
corresponding NiaAl-layer thickness. A clear oscillating 
behavior of the critical current density of more than 60 
single junctions versus dp is shown, indicating three dif- 
ferent - TT transitions. The amplitude of oscillations 
decays exponentially with a characteristic decay length 
of ^1=4.6 nm. It should be pointed out that the cor- 
responding oscillation period, given by 7r^2, is one mag- 
nitude smaller, namely ^2=0.45 nm. There is a rather 
good agreement of the theoretical fit after Eq. (^ with 
these two decay length's, see solid line in Fig. 21 The 
magnetically dead layer has been measured for different 
interfaces. Corresponding to this SIFS multilayer, we 
achieved a thickness of the dead layer around 5 nm to 
8 nm for each interface. So, we expected the first oscil- 
lation period at around 10 nm or even more. To find a 
theoretical explanation for the large difference between 
^1 and ^2 we start with linearized Eilenberger equations 
which are valid at the distances from SF interface larger 
than ^1 
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Here 6 is the angle between direction of the Fermi velocity 

and interface normal, t = vp/l is the scattering time. 

The solution of this equation has the form 



fix, 9) = Ci9) cxp {-^} , C// = ^r' + ^^2^'. 
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where £,eff is the effective decay length which is indepen- 
dent on 9 and C{9) is an integration constant. Substitu- 



tion of Eq. Q into Eq. |@J gives 
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Integration of Eq. (jHJ over angle 9 provides the equation 
for Ce// 
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2 .....X {In (1 + ^^0 ^) , In {i^H^) } , the solution of Eq. O 
reduces to Eqs. |(5J and 10), respectively. The results of 
numerical solution of Eq. Q are presented in Fig. [3 

There are steps on ^2 (Ch ) ^'^^ C2/C1 ^■^ ^h '^^~ 
pendencies (see Fig. [SJ accompanied by the minima on 
S,i^i£,H^) curves (see Fig. right inset). The ratio of 
^2/^1 falls very rapidly with increase of £,h- It follows 
from Fig. \S[ that the experimental value of this ratio 
^2/^1 ~ 0.1 can be achieved at I « lO^if. For the es- 
timated earlier parameters ^1 « 4.5 nm and ^2 ~ 0.5 nm 
from Fig.|Sl left inset and Fig.[Sl right inset we get nearly 
equal values for the electron mean free path I ^ 1.4^i ~ 
6, 3nm and I « 11^2 ~ 5 nm respectively. In the last esti- 
mation it was also supposed that l/£,o « 0.1, resulting in 
^0 ~ 50 nm and Fermi velocity vf = 1, 8-10^ m/s. With 
this value of vp and ^h ~ O.IZ w 0.5 nm we arrived at 
Eex ~ 1300 K « 0.11 eV. This combination of parame- 
ters is not unique. Starting from I « 20^h and l/^o ~ 0.3 
one can get I « 9.5 nm, ^0 ~ 30 nm, vp = 1.2-10^ m/s 
and Eex ~ 900 K fa 0,08 eV. These parameters are con- 
sistent with the previous experimental data integrated 
into Table H 

The discovered behavior of ^2 and ^1 is quite general 
and must be also observed in structures without ferro- 




FIG. 5: Dependence of the ratio ^2/^1 on inverse magnetic 
length l/S,H calculated for different ratios of l/i,o- Left inset: 
Inverse decay length Z/^2 vs inverse magnetic length 1/£,h for 
different ratios of Z/^o- Right inset: Inverse decay length Z/^i 
vs inverse magnetic length 1/£,h for different ratios of 1/S,q. 



magnetic ordering. An example is a normal filament of 
finite length, which is placed between superconducting 
banks and is biased by a dc supercurrent. It was shown 
|2C| . that the minigap induced to this filament from the 
S electrodes is not a monotonous function of phase dif- 
ference across the structure. This behavior could be also 
explained in terms of specific dependencies of ^2 and ^1 
upon electron mean free path in current biased systems. 
Summarising the presented results we conclude that 
utilization of the new ferromagnetic material, NisAl, as 
well as the wedge technique for its deposition permits for 
the first time the experimental demonstration as mush 
as six oscillations of the critical current as a function 



of thickness of ferromagnetic layer. High reproducibility 
of the junctions parameters, their scaling with the area, 
suppression of oscillation observed after only one change 
in operation-routing sequence, namely, replacement of Ar 
by Ne during the sputtering of NiaAl, clearly manifests 
that observed effect is due to magnetic ordering in NisAl 
film. The fact of this ordering has been also confirmed by 
independent examination of magnetization of the NiaAl 
films. It is important also to mention that NisAl is an in- 
termetallide. This metal is widely used and well studied 
before JL6|. It successfully combines the relatively small 
values of exchange integral with the transport properties 
close to that of strong pure ferromagnets. We believe 
that it will substitute the dilute ferromagnetic alloys in 
the SFS Josephson junction technology. The experimen- 
tal results are consistent with the theoretical predictions 
made in the frame of the Eilenberger equations. More- 
over, it was demonstrated that the intuitive knowledge 
about the relation between ^2 and S,, which is based on 
the dirty theories, has a very limited field of applications 
and can not be used for £^h > 51 or for E'e^T > 0.1. In 
particular, it was for first time recognized that the in- 
crease of Eex is not always accompanied by decrease of 
^1 sand there is some range of parameters when ^1 even 
may increase with Eex- The fact that one may combine 
reasonably large decay length with the smaller period of 
oscillations looks rather attractive for possible applica- 
tions of SFS Josephson junctions. 
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